Polaron Effects on Superexchange Interaction: 
Isotope Shifts of TJy, Tc, and T* in Layered Copper Oxides 
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A compact expression has been obtained for the superex- 
change coupling of magnetic ions via intermediate anions with 
regard to polaron effects at both magnetic ions and interme- 
diate anions. This expression is used to analyze the main 
features of the behavior of isotope shifts for temperatures of 
three types in layered cuprates: the Neel temperatures (Tn), 
critical temperatures of transitions to a superconducting state 
(Tc), and characteristic temperatures of the pseudogap in the 
normal state (T*). 
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Elucidating the nature of unusual isotope effect in 
copper-oxygen superconductors is one of the most impor- 
tant problems on the way to ascertaining the mechanism 
of the pairing of charge carriers in these compounds. It 
is known that the observation of an isotope shift of the 
superconducting transition temperature (Tc) was of cru- 
cial importance in ascertaining the phonon mechanism 
of pairing in conventional (low-temperature) supercon- 
ductors. The fact that an isotope effect exists in high- 
temperature superconductors built up of copperoxygen 
planes has long been beyond question; however, the rel- 
ative smallness of the coefficient (for example, a ~ 0.056 
for YBa2Cu4C>8 upon replacing 16 O with 18 O, instead of 
the standard value 0.5) and the specific features of its 
behavior in other compounds upon changing the number 
of holes in copper-oxygen planes do not fit the Bardeen- 
Cooper-Schrieffer (BCS) scenario. At the same time, the 
majority of authors of articles related to the isotope ef- 
fect (see the recent review [1]) correctly point out that, 
nevertheless, phonon modes, in some mysterious way, af- 
fect the superconducting transition temperature. In this 
context, we believe that the facts of observing an isotope 
effect for the characteristic temperature of the pseudogap 
state of underdoped cuprates (so-called pseudogap onset 
temperature T*) gain great importance. Thus, according 
to [2], the isotope exponent upon replacing 16 O with 18 O 
corresponding to T* is = 0.061. It was natural that this 
fact suggested a common origin of cct c and olt* [2]. 

We believe that another and even more important sim- 
ilarity in the dependence of the order parameters of the 
superconducting and pseudogap phases on the d - type 
wave vector (that is, cosq x — cosq y ) was explained in its 
time under the assumption that the transitions to both 
these phases are associated with short-range potentials 
[3]. It seems that superexchange interaction, screened 



Coulomb repulsion, and interaction of holes mediated by 
optical phonons are the most significant of these. The 
isotope shift of Tc and T* due to interaction mediated 
by optical phonons was discussed in a few works (see, 
for example [4,5]), and that for T* was considered in [6]. 
Below, we will focus our attention on polaron corrections 
to the superexchange coupling of copper spins (J) and 
will demonstrate that a number of features in the behav- 
ior of the isotope shift of Tc upon changing the num- 
ber of holes in the CuC>2 plane can be quite reasonably 
explained even within the framework of the purely su- 
perexchange mechanism of pairing. First, we emphasize 
the following important fact. Within the scenario [3], the 
superconducting transition temperature Tc ~ 2 J — G, 
whereas the characteristic temperature of the pseudogap 
phase T* ~ J + G [3,7]. Here, G is the parameter of 
the screened Coulomb interaction of holes on the nearest 
copper sites, which partially includes the correction due 
to interaction mediated by optical phonon modes. If it 
is granted that the isotope shift is associated with the 
phonon renormalization of G, the shifts of Tc and T* 
would be of opposite sign, which is in contradiction with 
the experimental results [2]. In this connection, it is be- 
lieved that the scenario of phonon renormalization of the 
parameter J is more reasonable to suggest as the source 
of positive isotope shifts of Tc and T* . 

The renormalization of J within the Hubbard model 
was investigated in detail by Kugel and Khomskii [8]. It 
is evident from the above estimates Tc ~ 2 J — G and 
T* ~ J + G that the small increase in J proportional 
to the phonon frequency found in this work gives the 
correct sign of isotope shifts but does not provide the re- 
quired magnitude of the shifts. This can be most simply 
demonstrated with the example of the isotope shift of the 
Neel temperature in related high-Tc compounds. As was 
already indicated in [1], the Kugel and Khomskii's cor- 
rection gives a correct sign of the isotope shift of the Neel 
temperature (T N ) in La 2 Cu0 4 upon replacing 16 O with 
ls O; however, even the most overrated estimates give a 
value that is six times lower than the experimental one. 

We believe that the main reason for the quantitative 
disagreement between the theory [8] and the experiment 
[1] is in the fact that the Hubbard model is not suit- 
able for the compounds of our interest. As was already 
stressed in [9] , that the energies of electron transfer from 
oxygen to a magnetic ion (A c ) in the majority of cop- 
per oxides are smaller than the energy of electron trans- 
fer from copper to copper (A a ). However, the Hubbard 
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model gives correct estimates only when A a > A c . 
The Hamiltonian in the form 

H = £ a a+a a + ^ £ c c+c<? + ^nfnj + 

+ J2 U cn c 1 n c l +J2t ac (a+c rJ + c+a a ). (1) 

is best suited to the description of superexchange interac- 
tion with taking into account explicitly the cascade hop- 
ping of electrons over oxygen. As applied to high- Tc su- 
perconductors, it is rather frequently named the Emery 
Hamiltonian. Here, t ac is the hopping integral between 
neighboring copper and oxygen sites, and U a and U c are 
the parameters of electron Coulomb repulsion. We will 
estimate the corrections to the superexchange parameter 
J due to polaron effects at copper (a) and oxygen (c) 
sites in the fashion of [8], supplementing Eq. (1) with 
the electron-phonon coupling operator 



(2) 



a,b,c 



Here, p q and V^ q are phonon annihilation and creation 
operators, and g, is the coupling parameter connected 
with the polaron stabilization energy (Ei) at site i by 
the equation Ei = g? jhiOi , where Wj are local vibrational 
frequencies. 

In the physical context, this calculation corresponds to 
a simplified Holstein model when the migrating charge 
locally interact with breathing modes, forming electron- 
vibrational states with dispersionless optical phonons. In 
this connection, it is pertinent to note that conduction in 
the compounds under consideration is just of the polaron 
type and is accomplished mainly via oxygen ion sites [10]. 

This calculation is most simply performed by the 
method of canonical transformations. The matrix of the 
unitary transformation of the initial Hamiltonian is found 
by excluding the odd terms with respect to hopping in- 
tegrals with an accuracy of to sixth-order perturbation 
theory. The calculation, whose mathematical details will 
be given in a more detailed article, gives the following 
result: 
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where A ac = e a — e c + U a — U c has the meaning of the 
energy of transfer from oxygen to copper, and the correc- 
tions proportional to E a hto a / A ac U^ and E a Tiuj a /U^ are 
not given because of their smallness for the compounds 
under consideration. Jo is the parameter of superex- 
change interaction of copper spins via the intermediate 
oxygen atom in the absence of phonons [9,11]. Note that 
the appearance of temperature factors in our equation is 



generally characteristic for the problems on transitions 
in transition metal compounds with the participation of 
quasilocal vibrations [12]. At the same time, one should 
keep in mind that polaron effects break down at T ~ u), 
and the concepts used here become inapplicable. 

The results of our calculation are given in the figure. 
The system of integral equations for the mean field pa- 
rameters corresponding to the transition to the pseudo- 
gap phase was solved self-consistently. We identify the 
pseudogap phase with the phase of sliding charge-density 
waves. This system was written in detail in [3] and is not 
given here. 

For the description of the superexchange coupling pa- 
rameter upon replacing some isotopes for other ones at 
Ucj/ksT, it is convenient to introduce parameters -feu 
and 70 by writing 
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It follows from Eq. (3) that, upon replacing 16 with 
18 0, 
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Substituting here (in electronvolts) Aqu-o = 1-5 [13] and 
standard values fiuj a — 0.05, E a = 0.4 [10,14], and us- 
ing the relationship AT/v/T/v ~ AJ/J characteristic for 
layered cuprates (see [1]), we find that the Neel temper- 
ature should decrease by 0.2% upon replacing 16 with 
18 in YBa2Cu3C>6.383- According to measurements in 
La2CuC>4 [15], the shift ~ 0.6%. If, however, it is as- 
sumed, following [1], that E a = 1.2 cV, our estimate will 
coincide with the experimental value. We hope that this 
explanation of the isotope shift of will stimulate fur- 
ther experimental investigations of this important prob- 
lem. Our estimated value 70 ~ — 0.014 is overrated. 
The value 70 ~ — 0.01 is better suited to comparison of 
the calculated with experiment, sec the figure. 

It is pertinent to note that, generally speaking, there is 
another possibility of changing J given by Eq. (4). This 
possibility is associated with the change in the distance 
between copper ions upon replacing some isotopes with 
other ones. It is known that the superexchange param- 
eters very strongly depend on the distance between the 
interacting ions. This mechanism explains well the in- 
crease in T c under the action of an external pressure on 
a high-T c crystal [16,17]. The question naturally arises 
in this case of what occurs with the lattice parameters 
upon changing some isotopes for other ones. Recent pre- 
cise measurements in YBa2Cu40s crystal showed [18] 
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that the lattice parameters a, 6, and c in the case of 
equal (in A) 3.8411(1), 3.8717(1), and 27.2372(8), respec- 
tively, whereas these are equal to 3.8408(1), 3.8718(1), 
and 27.2366(8), respectively, for 18 0, that is, these pa- 
rameters are somewhat smaller in the latter case. The 
positive isotope shift of the nuclear quadrupole resonance 
frequency of plane copper nuclei [18] is another important 
experimental fact, which indicates that interatomic cop- 
peroxygen distances are smaller in the case of 18 O. Based 
on these data, one may only conclude that the change of 
interatomic distances upon replacing some isotopes with 
other ones must lead to negative shifts of Tc and T* and 
will, probably, be relatively small. From the theoreti- 
cal point of view, this fact seems quite understandable, 
because changes in interatomic distances upon replacing 
some isotopes with other ones are due to the lattice an- 
harmonicity, and its effect is naturally of less importance 
than the effect of harmonic vibrations. 

0.5 . . . . . . 1 
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FIG. 1. The calculated exponent of the isotope shift of Tc 
upon replacing 16 O with 18 O as a function of the number of 
holes per one copper site. 

In the figure, the values of the coefficient ar c — 
-dln{T c )/d\n{M) for the replacement of 16 with 18 
are plotted as abscissas, and the numbers of holes per 
one copper site are plotted as ordinates. The symbols in 
the curve correspond to the points at which the system 
of self-consistent equations from [3] was solved. Because 
only the order of magnitude is known for the polaron 
energies E a and E c , the parameter 70 was normalized 
in such a way that ar c was equal to 0.1 at the optimal 
level of doping. The calculated behavior of ar c is not 
symmetric with respect to the point of optimal doping. 
The physical nature of this asymmetry is associated with 
strong competition between d-SC and irf-CDW phases in 
the underdoped state. On the left of this point, the iso- 
tope shift exponent increases with decreasing number of 
holes, approaching 0.5, whereas the value of oit c remains 
virtually constant as well as it does at 5/2 greater than 
0.16 (so-called overdoped regime). This kind of asym- 
metry (but without a step!) in the behavior of ar c as 



a function of the number of holes was found recently 
in measurements [19]. The authors of this work men- 
tioned already that, if only conventional interaction via 
the phonon field were responsible for the isotope effect 
and the unusual drop in the optimal doping region were 
related to a peak in the density of states, the curve would 
be approximately symmetric with respect to the point of 
optimal doping. At 8/2 larger than 0.16, the value of ar c 
would strongly increase; however, this was not found [19]. 

As to the isotope shift of Tc upon replacing the copper 
63 Cu isotope with 65 Cu or 66 Cu, the fact noted in [20] 
that the ratio aT c {Cu)/aT c (0) ~ 0.75 ±0.1 does not 
depend on the type of the compound and on the dop- 
ing level is naturally explained based on Eq. (3). The 
stabilization energy of a small-radius polaron (hole) at 
a copper site is higher than that at an oxygen site. The 
nearest environment of a hole at an oxygen site comprises 
positive copper ions, whereas the nearest environment of 
copper comprises negative oxygen ions. This fact is the 
reason for the difference between 70 and 7c«- 

The value of «t* calculated in this work in the re- 
gion 0.1 < 5/2 < 0.16 turned out to be approximately 
constant: «t« « 0.01. This is smaller than the value 
estimated in experiments (0.061) [2]; therefore, the effect 
of interaction via optical phonons on ay. discussed in 
[2] cannot be excluded. This is also corroborated by a 
number of experimental points in Fig. 2 from [19] on the 
right of the point of optimal doping. We hope to describe 
this problem in a more detailed work. 

Thus, the renormalization of the superexchange inter- 
action of copper spins due to polaron effects noted in this 
work explains the main regularities of the isotope shift 
of the superconducting transition temperature in layered 
cuprates both in the order of magnitude and in the sign 
and the character of the dependence on the number of 
holes. The starting equation for the renormalization of 
the superexchange parameter was verified using the iso- 
tope shifts of the Neel temperature of the parent com- 
pounds as an example. Our calculations were based on 
the scenario of competition between the superconducting 
phase and the charge-density-wave phase. Agreement 
between the calculations and experiment confirms this 
scenario. At the same time, our calculation predicts a 
rather sharp jump of the isotope exponent olt c on pass- 
ing through the point of optimal doping. This effect is 
relatively small; however, we believe that the experimen- 
tal observation of this effect would be of principal impor- 
tance. 
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